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Most evolution equations are partially integrable and, in order to explicitly inte- 
grate all possible cases, there exist several methods of complex analysis, but none 
is optimal. The theory of Nevanlinna and Wiman-Valiron on the growth of the 
meromorphic solutions gives predictions and bounds, but it is not constructive 
and restricted to meromorphic solutions. The Painleve approach via the a priori 
singularities of the solutions gives no bounds but it is often (not always) construc- 
tive. It seems that an adequate combination of the two methods could yield much 
more output in terms of explicit (i.e. closed form) analytic solutions. We review 
this question, mainly taking as an example the chaotic equation of Kuramoto and 
Sivashinsky vu'" + bu" + fiu' + u'^/2 + A = 0,u j^ 0, with (u, b, ^, A) constants. 



1. Introduction 

C ' Phenomena in continuous media are often governed by a partial differential 

^ , equation (PDE), e.g. in one space variable x and one time variable t 

X 



u: ^[\^">^-o- (" 



in which u and E are multidimensional, the integers to, n take a finite 
set of values. Our interest is the nonintegrable or even chaotic case, for 
which the powerful tools of Lax pairs, inverse spectral transform, etc ^ 
are inapplicable. The derivation of analytic results must then use other 
methods. Let us quote a few examples. 

(1) The one-dimensional cubic complex Ginzburg-Landau equation 
(CGL3) 

iAt + pA.^^ + q\A\^A - ijA = 0, pqlm{p/q) 7^ 0, (2) 

(and its complex conjugate, i.e. a total differential order four), in 
which p, q are complex constants and 7 a real constant, a generic 
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equation which describes many physical phenomena, such as the 
propagation of a signal in an optical fiber ^, spatiotemporal inter- 
mittency in spatially extended dissipative systems i^-iO'^s^ Yqy two 
coupled CGL3 equations, see analytic results in Ref. ^. 

(2) The Kuramoto and Sivashinsky (KS) equation, 

ift + VLpxxxx + h(pxxx + tJ.<Pxx + V^x =0, V ^Q, (3) 

in which v, 5, ji are real constants. This PDE is obeyed by the 
variable ip = arg^ of the above field A of CGL3 under some limit 
^^'^'', hence its name of phase turbulence equation. 

(3) The quintic complex Ginzburg-Landau equation (CGL5), 

lAt + pA^x + gl ApA + r| Al'^yl - i^iA = 0, pr Im{p/r) ^ 0, (4) 

in which p, q, r are complex constants and 7 a real constant. 

(4) The Swift-Hohenberg equation "^^'^^ 

iAt + bA^^^^ + pA^^ + q\A\^A + r\AfA - ijA = 0, br ^ 0, (5) 

in which b,p, q, r are complex constants and 7 a real constant. 

The autonomous nature of (1) (absence of any explicit dependence in x 
and t) allows the existence of travelling waves u = C/(^), solutions of the 
ordinary differential equation (ODE) 

u{x, t) = t/(e), ^^x-ct, E{U'^^\U^^-^\ ...,U',U)^0. (6) 

For the CGL3, KS, CGL5 and Swift-Hohenberg equations (with one 
exception, KS with 5^ = 16^i^), all the solitary wave solutions |Ap = 
/(^),(/3 — *&(C)i^ = X — ct, which are known hitherto are polynomials in 
tanhkS, (or cotanh, tan, cotan, which are the same in the complex plane), 
and such solutions are easy to find by taking advantage of the singularity 
structure of the PDE (see, e.g., the summer school lecture notes ^). 

Hence the natural questions: (i) Can other solitary waves u — f(x~ct) 
exist (in closed form)? (ii) If yes, please find them all, not just a few ones. 

The present paper introduces to the methods in principle able to answer 
both questions. They will mainly be exemplified with the KS equation (3). 

The paper is organized as follows. In section 2, we give a mathematical 
formulation of the problem. In section 3, we prove the inexistence of an 
analytic expression representing the general solution, and we compute the 
gap between the differential order N of the ODE (6) and the maximal 
number of integration constants in a singlevalued solution. In section 4, we 
give hints (not proofs) that some analytic result still has to be found. In 
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section 5, we review the consequences of the assumption of singlevaluedness 
for a solution of the ODE (6), and present an algorithm to implement them. 
In section 6, we present the consequences of the assumption of meromorphy 
for a solution of (6). The last section 7 states the open problems. 

2. Mathematical formulation of the problem 

The successive steps of the announced program are 

(1) To perform the traveling wave reduction from the PDE to an ODE. 
The KS PDE (3) depends on three fixed constants {v, b, y) (fixed 
means: which occur in the definition of the equation), the reduction 

ip{x, t) = c + u(^), S, = x-ct, (7) 

introduces in the ODE one more fixed constant A (the second con- 
stant c cancels out because of the Galilean invariance) 

vu'" + bu" + fiu' -\ \-A^O,i'^0, (8) 

and the general solution of (8), if it exists, depends on the four fixed 
constants (v, b, fx, A) and three movable constants (movable means: 
which depends on the initial data), which are the origin ^o of C ^-nd 
two other constants ci , C2 . 

(2) To count the number of constants which survive in the general so- 
lution of (8) when one requires singlevaluedness. 

(3) To find this largest singlevalued particular solution in closed form. 
Indeed, its representation as a series can be misleading, as shown 
by classical authors like Poincare and Painleve. 

3. Local separation of singlevaluedness and multivaluedness 

Because the ODE (6) is nonintegrablc, the number of integration constants 
present in any closed form solution is strictly smaller than the differential 
order of the ODE. This difference, an indicator of the amount of integra- 
bility of the ODE, can be precisely computed from a local analysis. 

Two local representations of the general solution of (6) exist. The first 
one, also the most well known, is useless for our purpose. This is the famous 
Taylor series near a regular point, whose existence, unicity, convergence, etc 
is stated by the existence theorem of Cauchy. The reason why it is useless 
is its inability to make a distinction between chaotic ODEs such as (8) and 
integrable ODEs such as u'" — 12mm' —1 = 0. 
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The second one, less known than the Taylor series of Cauchy, is a Lau- 
rent series (or more generally psi-series and/or Puiseux series) near a mov- 
able singularity xq. This one does provide the expected information. The 
technique to compute it is just the Painleve test (see Rcf. ^ for the basic 
vocabulary of this technique). Let us present it on the KS example (8). 

Looking for a local algebraic behaviour near a movable singularity xq 

U ^x^x„ Uqx'', Uo^O, X^ X -Xq, (9) 

one first balances the highest derivative and the nonlinearity, 

p-3 = 2p, p(p-l)(p-2)i.«o + Y=0, (10) 

a system easily solved as 

p^ -3, uq = 120i/. (11) 

The resulting convergent Laurent series, 

(o,_120^ 15b 15(16A^^-b^) 13(4M^-b> 

X"" X" 4 X 191.x 32x191.2 +^[X),(^^) 

lacks two of the three arbitrary constants. They appear in perturbation ^, 

u ^ «(") + e7i(i) + e2«(2) ^ ^ ^ ^ ^ (^3) 

in which e is not in the ODE (8). The linearized equation around u^"' 

-^+b^+f^^+-^"')-^''-0^ (14) 

dx"^ ax^ ax / 

has then the Fuchsian type near xq, with an indicial equation {q = —6 
denotes the singularity degree of the Ihs of (8)) 

limx-'-^i.dl+uoxnx''-' (15) 

= Hj - 3)(.? - 4)(j - 5) + 1201. = ,^{j + l){f - 13j + 60). (16) 

The resulting local representation of the general solution, 

u(a;o, ec_i,ec+, ec_) = 120i.x '^{Rcgular(x) 

+e[c-iX"^Regular(x) 

+ c+x(''+*^^/'Regular(x) 

+ c_x*""'^^^^'Regular(x)]+0(e')}, (17) 

in which "Regular" denotes converging series, depends on 4 arbitrary con- 
stants (xq, ec^i, ec^, sc_) but, as shown by Poincare, the contribution of 
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ec_i is the derivative of (12) with respect to xo, so c-i can be set to zero. 
The dense movable branching due to the irrational indices reflects ^^ the 
chaos, and to remove it one has to require £c+ = £c_ = 0, i.e. e = 0, making 
the analytic part of (17) to depend on the single arbitrary constant xq. 

The ODE (8) admits other Laurent series in the variable {u — y/—2A)^^, 
but they provide no additional information. 

The question is then to turn this local information into a global one, 
i.e. to find the closed form singlevalued expression depending on the maxi- 
mal number (here one) of movable constants. 

We will call unreachable any constant of integration which cannot par- 
ticipate to any closed form solution. The KS ODE (8) has two unreachable 
integration constants, the third one xq being irrelevant since it reflects the 
invariance of (8) under a translation of x. 

We will also call general analytic solution the closed form solution which 
depends on the maximal possible number of reachable integration constants, 
and our goal is precisely to exhibit a closed form expression for this general 
analytic solution, whose local representation is a Laurent series like (12). 

The above notions (irrelevant, unreachable) belong to an equation, not 
to a solution. Let us introduce another integer number, attached to a solu- 
tion, allowing one to measure its distance to the general analytic solution. 

The distance of a closed form solution to the general analytic solution 
is deflned as the number of constraints between the fixed constants and the 
reachable relevant constants. 

For the ODE (8), the fixed constants are i', b, jj,, A, the movable constant 
Xq is irrelevant, the movable constants ci = ec+, C2 ~ ec_ are unreachable, 
so the distance d is the number of constraints among the fixed constants. 

The closed form singlevalued solutions known to date are 



(1) one elliptic solution (distance d = I) ^'^'"^ 

2 , b/j, fi'^ 13^^ + vA 

in which p is the elliptic function of Weierstrass, 

p'^ = 4p3-52p-g3, (19) 

(2) six trigonometric solutions [d = 2) 16,13^ rational in e*^^, 

U - \2QVT^ - 15bT^ + ( — ^i - 30t.fc2 - ^^^' 



19 4x19j/ 

5,, 9 136^ 7ub k , k,^ , , , „, 

^2^^-32^719^ + 4^' ^--2^^-^-2(^-^o)m 
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the allowed values being listed in Table 1, 
(3) one rational solution (rf = 3), 

5 = 0, M = 0, ^ = : «= 120z/(e-^o)"^ 
which is a limit of all the previous solutions. 



(21) 



b'/iH 


uA/fi^ 


vk"^ / ^ 





-4950/19^, 450/19^ 


11/19, -1/19 


144/47 


-I8OO/473 


1/47 


256/73 


-4050/733 


1/73 


16 


-18, -8 


1, -1 



All those solutions admit the representation 



T) Log V' + constant, V = 60i^- 



15&- 



lb{l&^lu-b'^) d 



(22) 



d^3 ' ^""d^2 7gj^ d^' 

in which ip is an entire function. This linear operator I?, which captures 
the singularity structure, is called the singular part operator. 

The Laurent series (12) yields another information ^^. If its sum is 
elliptic, the sum of the residues of the poles inside a period parallelogram 
must vanish. Since the only poles of (8) are one triple pole, a necessary 
condition ^^ for the sum to be elliptic is to cancel the residue of (12), i.e. 
&^ = IQfiv. For this equation, the condition is also sufficient, see (18). 



4. Experimental and numerical evidence of missing 
solutions 

Experiments or computer simulations display regular patterns in the [x, t) 
plane (see ^3), some patterns being described by an analytic expression. For 
the other patterns, the guess is that there should exist matching analytic 
expressions. For the equation (3), one has observed a homoclinic wave ^^ 
Lp = /(^), ^ = a; — ct, while all known solutions are heteroclinic. 

The Laurent series (12) only provides a local knowledge of the general 
analytic solution. Rather than obtaining a global knowledge of the solution, 
which is the ultimate goal, it is easier to look at its singularities, by comput- 
ing the Fade approximants '^ of the Laurent series (12). Fade approximants 
are a powerful tool to study the singularities of the sum of a given Taylor 
series, and more generally to perform the summation of divergent series. 
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Given the first TV terms of a Taylor series near a; = 

N 



Sn^Yc.x^, (23) 



the Pade approximant [L, M] of the scries is the unique rational function 

T,m=0 b 



L,M]^ J^J-"^'''[ , bo^l, (24) 



obeying the condition 

Sn~[L,M]^0{x^+^), L + M^N. (25) 

The extension to Laurent series presents no difficulty. In particular, for L 
and M large enough, Pade approximants are exact on rational functions . 

The advantage of [L, M] over Sn (which has no poles) is to display the 
global structure of singularities of the series. 

From a thorough investigation of the singularities of the sum of the 
Laurent series (12) one concludes (this is not a proof): for generic values 
of {i/jb, ^, A), no multivaluedness is detected, no cuts are detected, and 
the singularities look arranged in a nearly doubly periodic pattern, the 
elementary cell being made of one triple pole and three simple zeroes. 

5. Consequences of singlevaluedness (Painleve) 

5.1. Classical results on first order autonomous equations 

The failure to detect any multivaluedness in the unknown general analytic 
solution by no means implies the singlevaluedness of this general analytic 
solution, because the Painleve test only generates necessary conditions, and 
the Pade approximants are a numerical investigation. It is however worth- 
while to examine in detail the consequences of an assumed singlevaluedness. 
Given the A^-th order autonomous algebraic ODE (6), any solution is 

«-./(e-eo), (26) 

in which ^o is movable. Provided the elimination of ^o between the equation 
(26) and its derivative is possible, one obtains the first order nonlinear ODE 

F{u,u')^0, (27) 

in which F is as unknown as /. 

However, /(^ — ^g) is now the general solution of (27), and there exist 
classical results on first order autonomous ODEs which are in addition 
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algebraic. Let us therefore assume from now on that the dependence of / 
on ^0 is algebraic (this is a sufficient condition for F to be algebraic). 

Let us summarize. Given the A'^-th order ODE (6) and its particular 
solution / Eq. (26), and assuming the dependence of / on ^o to be algebraic, 
one is able to derive a first order ODE (27) which is algebraic. 

Conversely, given an algebraic first order ODE F = Eq. (27), is it 
possible to go back to /? This question has been answered positively by 
Briot and Bouquet, Fuchs, Poincarc and put in final form by Painleve ^^ . 

Theorem 5.1. Given the algebraic first order ODE F ~ Eq. (27), if its 
general solution is singlevalued, then 

(1) Its general solution is an elliptic function, possibly degenerate, and 
its expression is known in closed form. 

(2) The genus of the algebraic curve (27) is one or zero. 

(3) There exist a positive integer m and (m + 1)^ complex constants 
aj,k, with flo^m 7^ 0, such that the polynomial F has the form 

F(m,u') = ^ ^ aj,ku^u' = 0, ao,„i 7^ 0. (28) 

Then, assuming / singlevalued with an algebraic dependence on ^Qj 

(1) It is equivalent to search for the solution / or for F. 

(2) The solution / can only be elliptic (i.e. rational in p and p'), or a 
rational function of e°^ with a constant, or a rational function of x. 

The explicit form (28) of F makes it much easier to look for F than /. 

5.2. Method to obtain the first order autonomous 
subequation 

The input data and assumptions are: 

(1) a A^-th order algebraic ODE (6), N >2, 

(2) a Laurent series representing its general analytic solution, 

(3) a first order algebraic ODE sharing its general solution with (6). 

Then, by the classical results of section 5.1, F has the form (28), and 
there exists an algorithm ^° yielding the solution / in the canonical form 

u = i?(p',p) = i?i(p) + p'i?2(p), (29) 
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in which _Ri, i?2 are two rational functions, with the possible degeneracies 
R{p',p)^R{e''i)^R{0, (30) 

in which R denotes rational functions. This algorithm is ^° : 

(1) Compute finitely many terms of the Laurent series, 

n = xnJ2 "^^^ + ^(X''+') J , X = e - Co- (31) 

(2) Choose a positive integer m and define the first order ODE 

„i [(m-fc)(p-l)/p] 

F{u,u') = J2 H aj^ku'u''' = 0, ao,my^O, (32) 

fc=0 j=0 

in which [z] denotes the integer part function. The upper bound on 
j implements the condition m{p — 1) < jp + k{p — 1), identically 
satisfied if p = — 1, that no term can be more singular than u'"*. 

(3) Require the Laurent series to satisfy the Briot and Bouquet ODE, 
i.e. require the identical vanishing of the Laurent series for the 
Ihs F{u.u') up to the order J 

F ^ x"<^-'^ ( E ^^X' + ^(^•'^') ) ' ^^ ■■ Pj = 0- (33) 

If it has no solution for a^jt, increase m and return to first step. 

(4) For every solution, integrate the first order autonomous ODE (32). 

The main step is to solve the set of equations (33), i.e. a linear, overdeter- 
mined system in the unknowns aj,fc. This is quite an easy task. 
An upper bound on m will be established in section 6. 

5.3. Results of the method on the KS equation 

The Laurent series of (8) is (12). In the second step, the smallest integer m 
allowing a triple pole (p = —3) in (32) is to = 3. With the normalization 
^0,3 = li the subequation contains ten coefficients, which are first deter- 
mined by the Cramer system of ten equations Fj = 0,j/ = : 6,8,9,12. 
The remaining overdetermined nonlinear system for [v, b, /i. A) contains as 
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greatest comnion divisor (gcd) b — IG/ii', which defines a first sohition 
— ~ ib, Us = u 



b V ( , b \ 9 f 2 15fe6 lOA^ ^ 



After division by this gcd, the remaining system for {ly, b, ii, A) admits four 
solutions (stopping the series at J = 16 is enough), namely the first three 
lines of Table 1, each solution defining a subequation, 

h n ,3^ 9 / 2 30m ,^30VV . (^a\ 

b^ 144 563 / 6 \3 9 4 „ 

— , U^ — U ;7, U -\ Us H u„ = 0, (37 

^iv 47 ' 144zy2 ' \ -iv ) mv " ' ^ ' 

b^ 256 456^ 



//i^ 73 ' 2048;/2 ' 



; b \ / , b \ 9/2 553 56 



(38) 

To integrate the subequations (34), (35)-(38), one first computes their 
genus '^, which is one for (34), and zero for (35)-(38). Therefore (34) has an 
elliptic general solution, listed above as (18). The general solution of the 
four others (35)~(38) is the third degree polynomial (20) in tanhfc(^ — ^o)/2- 

These four solutions, obtained for the minimal choice of the subequation 
degree m, constitute all the analytic results currently known on (8). 

6. Consequences of meromorphy (Nevanlinna) 

If the solution / is meromorphic, much can be said from the study of its 
growth at infinity (Nevanlinna theory). For the KS ODE, the meromorphy 
requires c+ = c_ = in (17), restricting the solution to the series (12). 
By direct application of the Nevanlinna theory, one can prove the 



'^For instance with the Maple command genus of the package algcurves ^^, which imple- 
ments an algorithm of Poincarc. 
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Theorem 6.1. ^ If a solution of (8) is meromorphic, then it is elliptic or 
degenerate of elliptic. Furthermore, 

(1) Elliptic solutions only exist if b^ = 16fiiy, and their order is three. 

(2) Exponential solutions have the necessary form P(tanhfc(^ — Co)); 
with k constant and P a polynomial of degree three. 

(3) The only rational solution is u ~ 120i^(C — Co) "^^ it exists for b =^ 
/i = ^ = 0. 

Consequently, the value m = 3 is an upper bound to the algorithm of 
section 5.2, which has therefore found all the meromorphic solutions of (8). 

7. Summary and open problems 

Let us represent the solutions of (8) by the following inclusions, 

elliptic C meromorphic C singlevalued C multivalued. (39) 

One has seen the various implications 

(1) (Singlevalued, algebraic dependence on xq) =^ elliptic (thm 5.1), 

(2) Meromorphic =^ elliptic (* using Nevanlinna theory), 

(3) Elliptic =^ (6^ = l&jiv) (residue theorem ^^), 

(4) (Elliptic or degenerate) => (order three) (^ using Nevanlinna the- 
ory) => (all such solutions in closed form ^°). 

The problem is open to find the general analytic solution in closed form 
for arbitrary [v, 6, jjl, A), which would be the sum of the Laurent series (12). 
Pade approximants and Painleve analysis find no multivaluedness nowhere. 

Two and only two possibilities remain about this general analytic solu- 
tion for generic values of [v, b, ii, A), 

(1) either it is multivalued, and strong efforts have then to be made to 
uncover this multivaluedness with both the Painleve test and the 
Pade approximants. This event is unlikely; 

(2) or it is singlevalued. In this case it cannot be elliptic, and the 
dependence on xq is necessarily transcendental. 

Solving this open problem would solve ipso facto many similar problems 
for nonintegrablc equations such as CGL3, CGL5 or Swift-Hohcnbcrg. 
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